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Recent STM experiments reveal niobium diselenide to support domains of striped (IQ) charge 
order side-by-side with its better-known triangular (3Q) phase, suggesting that small variations in 
local strain may induce a quantum phase transition between the two. We use a theoretical model of 
the charge order in NbSe 2 , based on a strong momentum- and orbital-dependent electron-phonon 
coupling, to study the effect of uniaxial strain. We find that as little as 0.1% anisotropic shift 
in phonon energies breaks the threefold symmetry in favour of a IQ state, in agreement with the 
experimental results. The altered symmetries change the transition into the ordered state from 
weakly-first-order in the 3Q case, to second order in the IQ regime. Modeling the pseudogap phase 
of NbSe 2 as the range of temperatures above the onset of long-range order in which phase coherence 
is destroyed by local phonon fluctuations, we find a shortening of the local ordering wavevector 
with increasing temperature, complementing recent X-ray diffraction observations within the low- 
temperature phase. 


The stability of different charge density wave (CDW) 
geometries is a topic of much current interest. In the 
simplest case a ‘IQ’ CDW selects a single preferential 
direction [T]. In two-dimensional materials, two perpen¬ 
dicular IQ CDWs can coexist, creating a 2Q pattern. 
In the layered high-temperature superconductors such a 
checkerboard state is believed to compete with IQ or¬ 
der HE]; recently the IQ state was found to dominate 
the 2Q in YBCO [6|. The layered hexagonal structure of 
niobium diselenide (2i/-NbSe2) allows instead for a 3Q 
triangular CDW with three superposed IQ patterns H- 
unj. Recent surface studies via scanning tunnelling mi¬ 
croscopy (STM) suggest a quantum phase transition be¬ 
tween 3Q and IQ states in NbSe 2 may be tuned by local 
strain, offering the possibility of studying the interplay 
between the two, as well as a direct comparison to the 
high-Tc case [TTl[T^ . 

It was recently shown that the minimal model repro¬ 
ducing the full range of experimental observations on the 
3Q ordered state in NbSe 2 includes a strong electron- 
phonon coupling, dependent on both the momenta and 
the orbital characters of the electronic states scattered 
between m- Motivated by the recent STM observa¬ 
tions, in this Letter we employ the same model in a free 
energy analysis to study the charge-ordered state upon 
application of uniaxial strain [T^]. We find that while 
the 3Q state, which respects the threefold lattice symme¬ 
try, is the lowest-energy configuration in the unstrained 
case, little uniaxial strain is required to break the sym¬ 
metry to IQ, agreeing with suggestions that lattice de¬ 
fects may locally stabilize the IQ state [12]. We find that 
the IQ phase transitions are second-order, whereas the 
3Q transitions are weakly-first-order. With increasing 
temperature we predict a shortening CDW wavevector 
in the locally-fluctuating, short-range-ordered, pseudo¬ 


gap regime above the CDW transition, analogous to the 
recently-observed non-monotonic evolution of the CDW 
wavevector within the ordered state |14j . 

Model-niobium diselenide consists of hexagonal layers 
of niobium atoms sandwiched between layers of selenium 
atoms, displaced so that they lie above and below half 
of the Nb interstitial locations. The crystal structure 
and bandstructure are depicted in Fig. Consecutive 
sandwich layers are displaced to have the complemen¬ 
tary half of the interstices occupied, giving two formula 
units per unit cell. 3Q charge order develops below 
Tcdw = 33.5 K. Recent STM measurements show sta¬ 
ble domains of IQ order on the surface of NbSe 2 , coex¬ 
isting with neighboring regions of 3Q order, and occur¬ 
ring in areas where the topmost niobium layer is slightly 
raised [T2]. This suggests that the IQ state may be sta¬ 
bilized by local strain, causing a quantum phase transi¬ 
tion between the IQ- and 3Q-ordered states [T2]. The 
CDW wavevector Qcdw = (1 ~ <5) fFM is found to be 
incommensurate with the lattice in both geometries, with 
SsQ « 0.014 and Siq « 0.143 (THHj. 

It has previously been shown that the CDW transition 
in NbSe 2 can be successfully described by a strong cou¬ 
pling between electrons if atomic displacements (/? 
through the interaction |I3j : 

Hint = X gk.k+q‘^q^k’^^k-rq: (1) 

where the electron-phonon coupling depends on 

the momenta of both the ingoing and outgoing electron 
states and the bands fi, v scattered between. The an¬ 
alytic expression for g has been shown to be applica¬ 
ble to a range of transition metal compounds includ¬ 
ing NbSe 2 , and accounts for the orbital characters of 
the electron bands with dispersion crossing the Fermi 
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FIG. 1. (a) The bandstructure of NbSe 2 modeled by the 

tight-binding fit described in Ref. 1131 Only the bands cross¬ 
ing Ey are shown, corresponding primarily to bonding and 
anti-bonding combinations of the two Niobinm or¬ 

bitals within a nnit cell. Only the lower band is involved 
in the CDW formation. (b) The Fermi surface consists 
of concentric barrel-shaped pockets aronnd the F and K 
points laiiaiii 116) . We ignore a small pocket around F. 
(c) The layered structure of NbSe 2 , with Nb in red and Se in 
blue, (d) The atomic displacements of the crystal structure 
caused by uniaxial strain within the plane of Nb atoms, in¬ 
dicated by the large blue arrows. Keeping the volume of the 
unit cell constant and increasing the horizontal interatomic 
distances to xo + e causes the diagonal interatomic distances 
to become approximately xo — e/2. 


level [lain]. It yields a stronger electron-phonon cou¬ 
pling in the lower band than the upper, explaining the 
order-of-magnitude difference in CDW gap size of the 
two bands in angle-resolved photoemission spectroscopy 
(ARPES) Haiiaiin]- We restrict attention to the lowest 
band from here on. This model reproduces the full host 
of available experimental observations on the 3Q state, 
with the overall strength of the electron-phonon coupling 
as the only free parameter (chosen to give the observed 
transition temperature Tcdw = 33.5 K) |13j . 

Following reference m we model the (isotropic) bare 
phonon dispersion Dg (|q|) by a Brillouin function with a 
maximum of 11.2 meV at the zone boundary. Since the 
longitudinal phonon is seen to soften in inelastic x-ray 
scattering [10], we restrict attention to g = q • g. Uni¬ 
axial strain is modeled as an anisotropic linear shift in 
bare phonon frequencies. As shown in Fig. [^d), stretch¬ 
ing the lattice in one direction causes a contraction in 
the perpendicular direction. We correspondingly define 
a strain parameter a through Di = Dg (1 -I- cr) along one 
of the Qcdw II TM directions, and ^ 2,3 = (l ~ f) in 
the 27r/3-rotated directions. 

Landau Theory-the free energy can be expanded in 
terms of the three order parameters ipj = {(p (qj)), where 
qj = (1 — (5) |rMj, and j labels the three inequivalent 
TM directions. At any given value for <5, the most general 


free energy to fourth order in pj is: 

3 

/3F [ipi,p2,P3] = - X! 1‘^j l^} “ 

i=i 

jmodS 

( 2 ) 

The phases of pj are chosen such that the third-order 
term is real and b is positive. Starting from the Hamil¬ 
tonian of Eq. Q, the coefficients a-d, and their depen¬ 
dences on q and (3 = 1/T, can be evaluated numeri¬ 
cally [3]. Strain enters via the Llj dependence of the 
quadratic coefficient. 

In the 3Q phase the order parameters obey ^ 1=^2 = 
= ip, while in each of the IQ phases only a single pj 
is nonzero. The free energy then takes the forms: 

I3F(q [‘Pj] = -a iqj,T, \pj\^ + c (q^, T) 

3 

/ 3 ^ 3 Q [^ 5 ] = - ^ a (q^-, T, Dj) \p\^ - b (q, T) 

+ 3c(q,T)(l + d(q,T))|(^|4, (3) 


where b-d do not depend on strain, and can be evaluated 
along any of the TM directions. Boundedness of the free 
energy requires c > 0, which is the case in all calculations 
described below. 

Starting from a disordered state with a < 0, mini¬ 
mization of Eq. @ yields a second-order phase transi¬ 
tion into either a 3Q or a IQ state whenever a becomes 
positive. If a stays negative, but c(l -|- d) < — a 3Q 
state forms via a first-order transition. Both scenarios 
require d > — 1, which is the case for all calculations con¬ 
sidered. For any nonzero 6, approaching a second-order 
transition into the 3Q state at a = 0 causes a diver¬ 
gence in b'^/a, which forces a first-order transition into 
the 3Q state instead. The two (would-be) transitions are 
typically close in temperature and strain, suggesting a 
weakly-first-order transition [3D]. We therefore only ex¬ 
pect either a second-order transition from disorder into a 
IQ state, or a first- or weakly-first-order transition from 
disorder into the 3Q state. 

Temperature Dependence-the quadratic coefficient in 
Eq. (§ describes the electronic susceptibility towards the 
formation of charge order, given by: 


a{qj,T, Llj) 
xiqj,T) 


-nj{qj)+x{qj,T), (4) 


k 


/(ek)-/(Ck-q) 

Ck — Ck-q — 


where / (^k) is the Fermi-Dirac distribution, and d —0+. 
The definition of the ‘generalized susceptibility’ x fol¬ 
lows Refs. nma and m and is expressed here within 
the Random Phase Approximation (RPA). The usual 
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(bare electronic) susceptibility, is given by x 

with g set to unity. Note that appear¬ 

ing in F 3 Q, is independent of strain. CDW order first 
sets in when a changes sign, meaning y has renormal¬ 
ized the phonon to zero energy at a nonzero wavevector 
q = Qcdw- At zero strain, Qcdw within RPA agrees 
with the experimentally-observed CDW wavevector m- 
To properly account for the system’s temperature de¬ 
pendence, it is necessary to include the entropy of the 
phonon field | 22 j . by introducing internal phonon lines 
into the RPA diagrams [23H25] . To fourth order in the 
coupling g, two such terms contribute to y, shown in 
Fig.i Their self-consistent inclusion forms the Mode- 
Mode Coupling Approximation (MMA), known to give 
more accurate temperature dependences than RPA in a 
range of transition metal compounds [53H1S]. We ne¬ 
glect one of the diagrams by appeal to Migdal’s theorem, 
since it contains only a vertex correction. The remain¬ 
ing diagram, constituting a self-energy correction to the 
electrons, is responsible for the observed suppression of 
the RPA transition temperature in NbSe 2 m- 

The range of temperatures in which RPA predicts or¬ 
der, while MMA predicts the order to be suppressed by 
the entropy of the phonon field, can be interpreted as a 
pseudogap regime m- In this regime the local ampli¬ 
tude of the displacements, y/{(p^(p), is nonzero, while the 
order parameter itself, ((^), averages to zero through fluc¬ 
tuations in phase [21H1S] . The result is fluctuating short 
range order without long-range phase coherence m- 
The CDW Wavevector-mxTueTicaWy evaluating the co¬ 
efficients in Eq. (§, within RPA, yields a weakly-first- 
order phase transition from a high-temperature disor¬ 
dered state to a 3Q CDW at Trpa = 416 K. Includ¬ 
ing phonon fluctuations via MMA suppresses the or¬ 
der, resulting in fluctuating regions of 3Q order with¬ 
out long-range coherence. At Tmma = 33.5 K we find 
a weakly-first-order transition in MMA, again to a 3Q 
state. This is interpreted as the onset of long-range 
phase-coherence EZ], as indicated in Fig.[| 

The CDW wavevector Qcdw coincides with the max¬ 
imum of the quadratic coefficient a (q^) in the free en¬ 
ergy, or equivalently of the generalized susceptibility y of 
Eq. Q. The shape of y is a convolution of the squared 
electron-phonon coupling with the bare electronic sus¬ 
ceptibility xo- The CDW wavevector is selected either 
by a divergence in xo due to Fermi surface nesting, by a 
strongly-peaked dependence of the electron-phonon cou¬ 
pling on electron momenta, or by a combination of the 
two; in NbSe 2 , Qcdw is selected predominantly by the 
wavevector dependence of g lilisi. 

While the MMA free energy expansion is strictly valid 
only close to the transition temperature, the susceptibil¬ 
ity can still be evaluated within the same approximation 
at temperatures above Tmma- Similarly the RPA expan¬ 
sion is valid only slightly above Trpa, but the susceptibil¬ 
ity can still be evaluated within the RPA-ordered regime. 



FIG. 2. The phase diagram of NbSe 2 as a function of uniax¬ 
ial strain and temperature. Cooling from high temperature, 
the disordered system develops fluctuating short-ranged order 
when the RPA predicts a transition, but long-range phase co¬ 
herence is suppressed by phonon fluctuations in the MMA. 
As the system cools further, long-ranged CDW order sets in 
at Tmma. The order changes from 3Q to IQ at an anisotropy 
in phonon energy of about 0.1%. Dashed black lines indi¬ 
cate second-order phase transitions, whereas solid black lines 
indicate weakly-first-order transitions. 

The peaks of x within these approximations provide an 
estimate for the CDW wavevector within the locally- 
fluctuating short-range-ordered region characterizing the 
pseudogap regime. As temperature increases, xo (Qj) 
broadens, and peaks in the temperature-independent 
electron-phonon coupling begin to dominate x- The re¬ 
sult is a decrease in the local | Qcdw I as temperature 
increases in the range Emma < Ecdw < Erpa, shown 
in Fig. A crossover between two distinct regimes, one 
close to Emma and the other at high temperatures, can 
clearly be seen. In the former regime, the bare suscepti¬ 
bility can modify the detailed momentum dependence of 
the electron-phonon coupling and cause small local ex¬ 
trema to appear in the full susceptibility, thus setting 
the precise value of Qcdw at low temperatures. Beyond 
the crossover, temperature has smoothed these local ex¬ 
trema, and the CDW wavevector is determined solely by 
the remaining momentum dependence of the electron- 
phonon coupling. This evolution complements recent X- 
ray diffraction observations within the ordered phase, of 
a non-monotonic evolution of Qcdw as temperature in¬ 
creases towards Ecdw from below M- 

Phase Diagram-tracking the coefficients of the free 
energy as functions of uniaxial strain and temperature 
yields the phase diagram of Fig. The range of tem¬ 
peratures between Erpa and Emma again corresponds 
to a pseudogap regime characterized by fluctuating lo¬ 
cal order [13 HZ]. Starting from unstrained conditions, it 
takes of order 0 . 1 % anisotropic change in the bare phonon 
energy to break the 3Q symmetry to IQ. To estimate 
the corresponding lattice deformations, we use the in¬ 
plane Griineisen parameter, — (a/Qo) dQo/da, with in¬ 
plane lattice parameter a. For several layered hexago- 
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FIG. 3. Left: the generalized susceptibility x along TM 
in the RPA at zero strain. From top to bottom, the black 
curve is at 33.5 K, blue 100 K, purple 303.5 K, and red 500 K. 
Right: the peak position of x, which dictates the CDW order¬ 
ing vector within the fluctuating, short-range-ordered regions, 
decreases with increasing temperature. Both the results of 
the MM A and RPA calculations are indicated. The solid and 
dashed lines are guides to the eye only. 


nal materials related to NbSe 2 , including M 0 S 2 , boron 
nitride, and graphite, the Griineisen parameter is of or¬ 
der unity [28l|29]. This suggests that a 0.1% anisotropic 
phonon energy shift also corresponds to changes of order 
0.1% in the inter-atomic distances. 

The transition temperature into the IQ phase increases 
with uniaxial strain, as the breaking of three-fold rota¬ 
tional symmetry stabilizes IQ order. The rates of in¬ 
crease at positive and negative strain differ because, as 
indicated in Fig.[^d), the atomic displacements stabiliz¬ 
ing the IQ phase are proportional to tr at positive strain, 
when the CDW order develops along the strain direc¬ 
tion, but proportional to —(t/2 for negative strain, with 
CDW order at an angle of 27r/3 to the direction of ap¬ 
plied strain. We found the 2Q state, with two CDWs at 
±27r/3 to the strain direction, to have higher free energy 
than either the IQ or the 3Q state within the range of 
parameters relevant to NbSe 2 [50] . 

At Tcdw = Tmma and zero strain, the system orders 
into a 3Q geometry via a weakly-first-order transition. 
Within the MMA, phonon fluctuations renormalize the 
generalized susceptibility, yielding a decreased range of 
stability for the 3Q state as compared to the RPA. Ap¬ 
plication of just under 0.1% strain then stabilizes the IQ 
CDW, in agreement with STM findings of a IQ phase 
in regions with an upper-bound uniaxial strain of about 
0.45% [12]. We predict this IQ order to also be stabilized 
in the bulk upon application of uniaxial pressure to the 
entire sample. 

Conclusions-We compared the stability of proximate 
IQ- and 3Q-ordered CDW phases in niobium diselenide 
using a Landau free energy expansion, based on a tight- 
binding model that was previously shown to capture all 
of the various known properties of the zero-strain, 3Q- 
ordered phase m- By comparing the model’s phase dia¬ 
gram within the RPA and MMA, it is found that phonon 
fluctuations, present in the latter, suppress long-range 
order and lead to the formation of a pseudogap regime 
between Trpa and Tmma in which the order parameter 


has a nonzero local magnitude but no long-range phase 
coherence. Both the pseudogap and CDW phases are 
shown to have a 3Q geometry, but it takes only of or¬ 
der 0.1% uniaxial strain to induce a IQ geometry, in 
agreement with suggestions that the IQ state seen in 
STM is stabilized by surface-layer distortions causing lo¬ 
cal strain [T5] . We propose that the bulk quantum phase 
transition between the different ordering geometries can 
be accessed by the application of comparable amounts of 
bulk uniaxial strain. 

The free energy expression shows that while, in prin¬ 
ciple, the 3Q phase transition can be either first- or 
second-order, any would-be second-order transition is 
pre-empted by a weakly-first-order transition. We found 
that the 3Q phase at low strain is of this weakly-first- 
order form, in agreement with previous experimental and 
theoretical studies [23 [5nU55] , while the transition from 
disorder to the strained IQ phase is second order. 

The CDW wavevector Qcdw is dictated by the loca¬ 
tion of the maximum in the momentum-dependent gen¬ 
eralized susceptibility x, implying that within the short- 
range-ordered domains of the pseudogap, Qcdw can be 
tracked by the peak in x between Tmma and Trpa. As 
temperature increases within the pseudogap, the momen¬ 
tum dependence of the electron-phonon conpling increas¬ 
ingly dominates the bare electronic susceptibility in de¬ 
termining the shape of %, resulting in a decrease of the 
CDW wavevector within the locally-fluctuating regions. 
This complements recent experimental work in which a 
non-monotonic evolution of the CDW wavevector was 
fonnd within the ordered state itself El- 

The situation in NbSe 2 recalls that in the layered high- 
temperature superconductors, which exhibit both charge 
order [55H55] - believed to be close to the phase boundary 
between the IQ and 2Q geometries supported by their 
square lattices - and pseudogaps characterized by 
some form of fluctuating local order EHSmO]. Contrary 
to the high-Tc case, the underlying physics giving rise 
to these features in NbSe 2 is well understood [T5]. The 
interplay of the CDW geometries with external param¬ 
eters such as applied uniaxial stress may now be used 
as a test case from which to consider the geometries of 
charge-ordered states more generally, and may thus shed 
some light on the case of the high-Tc superconductors. 
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